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1 Introduction 

Consider the Galton-Watson branching process with immigration, {y{k) : k = 1,2, ■ ■ ■ }, defined by 

y{k-l) 

y{k)= Yl C{k,j)+r]ik), k>l, y(0)=0, (1.1) 
i=i 

where {S,{k,j) : k,j = 1, 2, • • • } and {7?(A;) : k = 1,2,- ••} are two independent families of i.i.d. 
random variables taking values in N := {0, 1,2,---}. The distribution of i{k,j) is called the 
offspring distribution and the distribution of i]{k) is called the immigration distribution. Let g{-) 
and /i(-) be the generating functions of i{k,j) and rj{k), respectively. It is easy to see that {y(fe)} 
is a discrete-time Markov chain with values in N and one-step transition matrix P{i,j) given by 

oo 

J2 P{h 3)s' = 9{sYh{s), i G N, < s < 1. (1.2) 

i=o 

For simplicity, we also call {y(A;)} a GWI-process with parameters (g, h). Assume that the offspring 
mean m := g'{l) is finite. The cases m > 1, m = 1 and m < 1 are referred to respectively as 
supercritical, critical, and subcritical. A sequence of GWI-processes {yn{-)} with {gn,hn) is said 
to be nearly critical if m„ := (7^(1) converges to 1 as n tends to oo. 

The estimation problem for the offspring and immigration parameters in the GWI-process has 
been extensively studied; see Heyde and Seneta [SI [U], Wei and Winnicki and the references 
therein. It is well known that the conditional least squares estimators (CLSE), first obtained by 



Klimko and Nelson [13] , can be used to estimate the offspring mean on tlie basis of the observing 
information on {y{k)}; see also [8] and [25] for other closely related estimators. In the non-critical 
case, the CLSE of the offspring mean is consistent and asymptotically normal (see [H], [241 125j). 
However, it was shown by [22] and [23] that in the critical or nearly critical case the CLSE is 
not asymptotically normal. In fact, when the process is nearly critical and the offspring variance 
tends to a positive real number, Sriram [22] gave the weak convergence of GWI-processes to the 
branching diffusion with immigration. As a result, the above CLSE of the offspring mean has the 
asymptotic distribution which is expressed in terms of the limit process and the normalizing factor 
is n. Motivated by the similar statistical application, Ispany et al. [TT] have recently obtained a 
fluctuation limit theorem for the nearly critical processes where the offspring variances tend to 
0. Such limit is a time-inhomogeneous Ornstein-Uhlenbeck (OU) processes driven by a Wiener 
process. As a consequence, they proved the asymptotic normality of CLSE of the offspring mean 
with normalizing factor n^/^. Obviously, the asymptotic behavior of the CLSE in the critical or 
nearly critical case is closely related to the limit theorems of the GWI-processes. 

The main objective of this paper is to give a general fluctuation limit theorem and its applications 
for processes that allow the offspring and immigration distributions to have infinite variances. 
Fluctuation limits for branching models with immigration have been investigated by Dawson and 
Li [1], Ispany et al. [IT], Li [18] and Li and Ma [20] : see also Dawson et al. [3] for the type 
of limits in the measure-valued setting. In the present paper, we shall consider a sequence of 
nearly critical GWI-processes {yn{-)} with {gn,hn) satisfying a set of conditions similar to that 
of [18]. Let us define the sequence Yn{t) = y„([nt]) and consider the rescaled centralized process 
Zni') = c~^(y„,(-) — with certain sequence of positive constants c„. It turns out that 

Zni') converges to a time-inhomogeneous OU type process driven by a spectrally positive Levy 
process (Theorems 12. 2p . Based on this fluctuation limit, we show that non-degenerate limit laws 
still exist for the above CLSE estimates of means (Theorem 13. 2p . Of special interest is the case 
when the offspring and immigration distributions belong to the domain of attraction of a stable 
law with exponent a (1 < a < 2). For simplicity, suppose that 

5n(s) =s + -(l-s)" and K{s) = s + w{l - s)" , 
n 

where < 7, < 1/a. Note that for 1 < a < 2^ has infinite variance but its heavy-tailed 
effect weakens as n — > 00; for a = 2, the offspring variance is 27/n and tends to 0. Then Zn{-) 
with Cn = n^/" converges to a OU type process driven by a a-stable process (Corollary 12. 3p . As 
a consequence, the CLSE of the offspring mean is asymptotic to a a-stable distribution and the 
normalizing factor is (Corollary [XT]). As mentioned above, the estimation for the offspring 

mean in GWI-process have been systematically studied by [241 [25], [22] and [llj . provided that 
the offspring variances are finite. Our results can be regarded as an attempt in the case when the 
above assumption fail to hold. 

Another interesting case, related to our limit theorem, is that the offspring variances are finite 
and tend to 0, but the offspring distributions do not satisfy the Lindeberg conditions required in 
|llj . Then the resulting fluctuation limit Z(-) is a OU type process with positive jumps instead 
of OU diffusion (Corollarv I2.2p . In this case, it is also possible to consider the CLSE estimates 
for the offspring and immigration variances. We show that the CLSE of the offspring variance 
is consistent and its asymptotic distribution (with normalizing factor n) is expressed in terms of 
the jumps of Z{-), while the CLSE of the immigration variance is not consistent (Theorme 13. Sp . 
However, if we return to the case of [11] (see also Example 12. ip . the above asymptotic distribution 
is degenerate to and the above immigration variance estimator becomes consistent (Remark 
13. ip . Hence, in this case, by adding certain conditions on fourth moments we further prove that 
these estimators of the offspring and immigration variances are asymptotically normal with the 
normalizing factors ri?^"^ and n^^^, respectively (Theorem 13. 4p . This result also contrasts with 
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the critical-mean and positive- variance case of Winnicki [26] . in which the CLSE of the offspring 
variance is not asymptotically normal, although it has another limit law with the normalizing 
factor n}!'^ . 

The remainder of this paper is organized as follows. The main limit theorems and some exam- 
ples will be given in section 2. In Section 3 we obtain some asymptotic estimates for the statistics 
of the GWI-process, as applications of our limit theorems. Section 4 is devoted to the proofs of 
Theorem E^ljlS^] and Theorem 01331 

Notation. Let M+ = [0, oo). For x G M, set x{x) = (1 A x) V (-1). " ^ " and " ^ " denote 
the convergence of random variables in probability and convergence in distribution, respectively. 
We also make the convention that J^* = ~ // = /(.^ and = /j.^, for r < t G M. 

2 Limit theorems and examples 

Let us consider a sequence of GWI-processes yn(') with parameters ((7„, A realization of ?/„(•) 
is defined by 

ynik)= Uk,j)+Vn{k), k>l, yn(0) = 0, (2.1) 

i=i 

where {Cnik,j)} and {r]nik)} are given as in (jl.ll) . but depend on the index n. Also, gn and hn 
are the generating functions of £,nik,j) and r]n{k). Now introduce the sequence 

Yn{t) := yn{[nt\), t > 0, 

where [nt] denotes the integer-part of nt, and Y^{t) := Ylt^=i Vn{k)- We first prove a limit theorem 
for the sequence (l^(-), 1^(-)). Such theorem is the modification of Theorem 2.1 in [19]. Let {hn} 
be a sequence of positive numbers such that — *■ oo as n — > oo. For < A < 6„, set 

Rn{\) = nbn[{l - X/bn) " 5n(l " A/5„)] (2.2) 

and 

FniX) = n[l - Kil - X/bn)]. (2.3) 
Consider the following set of conditions: 

(A) The sequence {Rn} is uniformly Lipschitz on each bounded interval and converges to a 
continuous function as n — > oo; 

(B) The sequence {Fn} converges to a continuous function as n ^ oo. 
Lemma 2.1 Under condition (A), the limit function R of {Rn} has representation 

/•oo 

R{X) = cX-9X^- / (e~^" -l + Xu) Ai{du), (2.4) 
Jo 

where c € M, 9 > 0, and Ai{du) is a a-finite measure on (0,oo) with J^{u Au'^)Ai{du) < oo. 
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Lemma 2.2 Under condition (B), the limit function F of {Fn} has representation 

POO 

F{X) = d\+ (1 - e-^")A2((in), (2.5) 
Jo 

where d> 0, and K2{dz) is a a-finite measure on (0, oo) with /q°°(1 A u) A2{du) < oo. 

Theorem 2.1 Suppose that (A) and (B) are satisfied. Then (Yn{-) /bn,Y^{-)/bn) converges in 
distribution on D{[0, oo), M^) to a two-dimensional non-negative Markov process [Y{-), Y'{-)^ with 
initial value (0,0) and transition semigroup {Pt)t>o given by 

j e-<^'">Pj(x,dn) =exp| -xi^j(zi)-X2Z2- /" i^(V's(^i)+^2)ds|, (2.6) 

where x = (xi, X2) G z = {zi, Z2) £ and V't (-^i) is the unique solution of 

^{z,) = RiMzi)), Mzi)=zi. (2.7) 

Remark 2.1 (i) The process Y{-) is a conservative continuous state branching process with im- 
migration (CBI-process) and Y'{-) is the immigration part of Y{-). See Kawazu and Watanabe 
[13j for a complete characterization of the class of CBI-processes. Furthermore, {Y{-),Y'{-)) is a 
special case of two-dimensional CBI-processes; see |21] and the references therein. 

(ii) Lemma 12.1112.41 are closely related to the Levy-Khinchin type representations of some class 
of continuous functions. Li [1^ used a simple method based on Bernstein polynomials to prove 
Lemma l2.ll and 12.21 Here, inspired by Venttsel' [23], these lemmas can be proved in a different 
but more intuitive way. From the proof of Lemma |2.3^ we can see that either (A) or (D2) implies 
the convergence of the sums of the triangular array of i.i.d. variables {^lilMLJ. ; k,j = 1,2- ■ ■} or 
)^ inikj)-i . j^^j _ i 2---^. See Grimvall [7] for a similar consideration. 

Corollary 2.1 (IJ^jf) Let L{x) and L*{x) be positive functions slowly varying at 00 such that 
L{x) ~ L*{x) as X — > 00. Consider a sequence of GWI-processes with {gn, hn) given by 

gn{s) EE s + 7(1-s)"l(^), 

hn{s) = l-ti7(l-s)"-iL*(^), 

where I < a < 2, 7 > 0, tx? > 0. Let bn be the sequence satisfying 

bn ~ [nL(6„)]V(--i) (^ [nL*(6„)]i/("-i)). (2.8) 

Then {Yn{-) /bn,Y^{-) /bn) converges in distribution on L»([0, cx)), M^,) to {Y{-),Y'{-)) defined by 

dY{t) = ^Y{t-)dX{t) + dY'{t), y(0) = y (0) = 0, (2.9) 

where X{t) is a spectrally positive a- stable Levy process with Laplace exponent R{X) = — 7A", 
and Y'{t) is a {a — 1)- stable subordinator with Laplace exponent F{X) = wX"~^, independent of 
X. 

Remark 2.2 Sometimes the above process can be regarded as the branching process conditioned 
on not being extinct in the distant future, or Q-process; see Lambert |16] . The pathwise uniqueness 
for the type of SDE ()2.9p has recently been proved by Fu and Li [6]. 
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Proof of Corollary \2.1\ Without loss of generality, consider Rn{X) by (12. 2|) on A G [0, 1]. It is 
easy to see that lim^^oo -Rn(A) = -7A". For A > 0, \R'JX)\ = a7A"-^[nL^(6„/A)/6^-^], where 
L'^{x) is a slowly varying functions such that L'^ix) ~ L{x), as x — > 00 (cf. j2l Theorem 1.8.2]). 
By the representation theorem of L'^{-) ([21 Theorem 1.3.1]), 

r /-fen/A 

L%bn/\)/L\hn) = {q{hn/\)/q{bn)} exp / e{u)du/u 

where q{x) q £ {0, 00), e(x) — > as x ^ 00. Fix < e < g A (a — 1). There exists xq > such 
that q — e < q{x) < q + e and \e{x)\ < e, if x > xq. Note that 6„/A > 6„ for < A < 1 and choose 
sufficiently large n, we have 

A°-i[L^(6„/A)/LO(6„)] < i±£A°-i-^ 

q-e 

Then, by the above inequality and ()2.8p . sup„ ]i?^(A)] is bounded in A G (0, 1]. Note that i?^(0) = 
and thus (A) holds. Also, it is not hard to see that lim„^oo 

(A) = roA°"^ By Theorem ED the 
limit process {¥{■),¥'{■)) is defined by I^B) and ([221) with R{X) = -7A° and F(A) = tzjA"^^ By 
[6], (Y{-),Y'{-)) is the unique solution of the above stochastic equation system. □ 

Now we turn to study the fluctuation limit for the sequence Yn{-). Assume that m„ = 5^(1) 
and LOn = /in(l) finite. Let {c„} be a sequence of positive numbers. Set 

G„(A) = n^iil - m„A/c„) - 5„(1 - A/c„)] (2.10) 

and 

Hn{X) = n[(l - tOnX/Cn) " K{1 - A/c„)], (2.11) 
for < A < c„. We will need the following conditions: 

(C) n/c„ — > oo and n/c^ — > 70 as n ^ 00, for some 70 > 0; 
(Dl) n{mn — 1) ^ a as ?i — > 00, for some a G M; 

(D2) The sequence is uniformly Lipschitz on each bounded interval and converges to a 

continuous function as n — > 00; 

(El) a;„ — > w as n ^ 00, for some a; > 0; 

(E2) The sequence {Hn} is uniformly Lipschitz on each bounded interval and converges to a 
continuous function as n ^ 00. 

Lemma 2.3 Under conditions (C) and (Dl,2), the limit function G of has representation 

/•oo 

G(A) = /3iA- CTiA^ - / (e"^" - 1 + An)/i(dn), (2.12) 
Jo 

where /3i G M, (Ji > and 2ai > 070, and fi{du) is a a-finite measure on (0,oo) with J^{u A 
u'^)fi{du) < CO. 
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Lemma 2.4 Under conditions (C) and (El, 2), the limit function H of {Hn} has representation 

/•oo 

if(A) = /JsA - daA^ - / (e-^" - 1 + An)z^(dn), (2.13) 
Jo 

where G I^, c"2 ^ and 2a2 + w7o > <^^70i and i^(du) is a a-finite measure on (0, oo) with 
J^{u A v?')v{du) < oo. 

Let <j){s) = u) Jq* e'^^du and let q{s) = {2ai — a7o)</'(s) + 2o"2 + a;(l — u;)7o- By the above represen- 
tations, g{s) > for s > 0. We actually obtain a set of parameters (/3i, /?2, £>(•))</'(■); A*) ^) which 
will be used to characterizes our limit processes. Let Zn{-) be defined by 

Z„«) = (2.14) 
Our main result of this paper is the following fluctuation limit theorem. 

Theorem 2.2 Suppose that conditions (C), (Dl,2) and (El, 2) are satisfied. Let B{t) be a one- 
dimensional Brownian motion, No{ds,du) be a Poisson random measure on (0, oo) x M_|_ with 
intensity dsv(du) and Ni(ds,du,dC) be a Poisson random measure on (0, oo) x M+ x (0,oo) with 
intensity dsfj,{du)dC. Suppose that B, Nq and Ni are independent of each other. Then Zn{-) 
converges in distribution on D([0,oo),M) to a time-inhomogeneous OU type process Z(-), which 
can be constructed as the unique solution of the following stochastic equation 



Z{t) = / {l32 + l3i(t>{s) + aZ{s))ds+ y^dB(s) 
Jo Jo 

ft r ^ ft f f<P{s) 

+ uNo{ds,du)+ / / uNi{ds,du,dC), (2.15) 



wiiere iVo((is, du) = NQ{ds, du) — dsv(du) and Ni{ds, du, dQ) = Ni{ds, du, d() — dsfi{du)d(. 

Remark 2.3 The conditions of Theorem 12.21 imply that Yn(t)/n converges weakly to the deter- 
ministic function cj){t) = u> /q e'^^ds. In fact, consider i?„ in (12. 2p and F„ in (12. 3p with 6„ = n. Note 
that = Gn{cnX/n) + n{mn — 1)A. Then by conditions (C), (Dl,2) and Lemma [TTl i2^(A) 

is uniformly bounded in each bounced interval and lim^^oo -Rn(A) = a\. In a similar way, we also 
have lim„^oo -^ri(A) = oo\. The above weak convergence result follows from Theorem 12.11 

Corollary 2.2 Let {y„,(fc)} be defined as in the beginning of this section. In addition to conditions 
(Dl) and (El), we assume that 7r„, = var,^„(l,l) < oo, r„ = varr/„,(l) < oo, and the following 
conditions hold: 

(a.l) The sequence fLn{-) = nE (^„(1, 1) - "in)^l{(5„(i,i)_m„)/v^ e ■ } converges weakly to a finite 
measure denoted by /i(-), as n — > oo; 



(a.2) The sequence P„(-) = E^(7?„(l) - ujn)^l^^^^^-^)_^^y^ g.^ 
sure denoted by i>{-), as n ^ oo. 



converges weakly to a finite mea- 



Let Zn{-) be defined by i2.14\) with Cn = \/n. Then Zn{-) converges in distribution on D{[0, oo),] 
to a OU type process Z{-) whose Levy measure has finite second moment, i.e. 

Z{t) = a f Z{s)ds+ I ^/^dB{s)+ I I uNo{ds,du) 



+ / / / uNi{ds,du,dC), (2.16) 
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where g{s) = t'({0}) + /l({0})(/>(s) and similarly cl){s) = lo Jq e°''^du. B, Nq and Ni are defined as in 
Ii2.15\) . hut the corresponding intensities of Nq and Ni are given hy dsv{du) and ds^{du)dC, with 
v{du) = ii~2l|u>o}^('^^) 3.iid fj,{du) = u~'^l^uyQjfl{du). 

Proof. By Theorem 1 2. 2 1 it suffices to check (D2) and (E2) are satisfied. First it follows from (a. 1,2) 
that /i(-) and D{-) are supported by [0,oo). Consider G„,(A) in (|2.10p with c„ = ^/n. Without loss 
of generality, we restrict ourselves to A E [0, 1]. Making a Taylor expansion of gn about 1, we have 
G„{\) = -n J^{1 - s)g'^{l - s\/^)\^ds. Note that 

+ nm„(m„ - 1)E[(1 - sA/^/^)^"^^'^)-'""] 
+ n(2m„ - l)E[(e„(l, 1) - m„)(l - sA/V^)«"(i'i)-™" 

Fix s, A € (0, 1] and choose sufficiently large n such that < s\/ ^/n < 1/2. We have that 
^nv^mfi-Vv^)^^^^^^ _ T-nsA-^^^^) < /(e^"^^ | n| / V^) V (e(e _ 1) (rf^) , 



|l - (1 - sA/^A^)«"(^'i)-"^" I < 2e2|e„,(l, 1) - mn\/V^. 

Note that fin{-) is supported by {{k — mn)/^/n : k = 0, 1,---}. By conditions (Dl), (El) and 
(a.l), it is easy to see that ng"{l — s\/y/n) — > a + /[q^^-) e~^^^fi{du). Since ng"{l — sX/^/n) is 
bounded, (D2) holds and lim„_,oo G„(A) = -(/i({0}) + a)X^/2 - //ooo)(e~^" " ^ + >^u)/u^Kd' 



u). 



It follows in a similar way that (E2) also holds and lim„^oo -f^n(A) = — (^'({0}) + oj — u;)A /2 

f(0,oo)( 



-L^.{e-^''-l + Xu)/u^i>{du). □ 



Example 2.1 Theorem 2.2]) Assume that (Dl), (El) and the following conditions hold: 

(b.l) niTn TT and r„ ^ r as n ^ cxd for some vr > and r > 0, 

(b.2) nE (^„(1, 1) - mnfl{i^„(iA)-m„\>^e ^ as n ^ oo for all e > 0, 

(b.3) E (r/„,(l) - a;„)2l||^^(i)„^,^|>^^} ^ as n ^ cx) for all e > 0. 

In this case, we see that conditions (a.l, 2) are satisfied with fl{du) = 7r(5o(dn) and i){du) = r6o{du) 
{5x{du) denote the dirac measure at u = x). Then we still have the above limit theorem, and the 
fluctuation limit process Z(-) is given by 

dZ{t) = aZ{t)dt + y^dB{t), Z(0) = 0, (2.17) 

where g^t) = r + Tr<j){t), and B{t) is a one-dimensional Brownian motion. 

Example 2.2 Suppose that P(^„,(l, 1) = [^/n]) = l/n^ and P(^„(l,l) = 1) = 1 - l/n^, while 
P(??n(l) = iV^]) = '^/n and P(r/„(1) = 1) = 1 — 1/n. We see that n(m„ — 1) ^ 0, u;„ ^ 1 and 
(a.l, 2) are satisfied with jl{du) = D{du) = 6i{du). Another example is as follows. Suppose that 
fl{du) is any non-degenerate finite measure on (0, oo). For large enough n, let iJ,{du) = u~'^jl{du), 
Unidu) = /i((l/ni/^,cx3))"^l|„>„i/4}/i(dn), and 

gnis)=Pn r e-^^^'-^^fin{du) + {l-pn), where pn = r 
Jo /Li((l/nV4^cx))) 
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Let ^„(1,1) have the distribution corresponding to (/«(•)■ Note that f^^-^-^^^^l/4yUfi{du)/^/n — > 

and /i((l/n"'^/^, oo))/n 0. Then it is not hard to see that n{mn — 1) ^ and condition (a.l) is 
fulfihed with ji{du). r]n{l) can be constructed in a similar way. 



Corollary 2.3 Let L{x) and L*{x) he positive functions slowly varying at oo such that L(x) 
L*{x) as X — > oo. Consider a sequence of GWI-processes with (gn, hn) given by 

gn{s) = (1 - mn) + mnS + -(1 - 



n \1 — s 

hn{s) = (1 - UJn) + (^nS + w{l - s)"" L* f-^— , , 

\ \ — sJ 

where 1 < a < 2, 7 > 0, > 0. m„ and u;„ satisfy conditions (Dl ) and (El). Let Zn{-) he defined 
by Ii2.14\) with Cn satisfying 

Cn ~ [nL(c„)]^/" (~ [nL*(c„)]V-). (2.18) 
Then Zn{-) converges in distribution on D{[0, 00), M) to a OU type process Z{-) defined by 



dZ{t) = aZ{t)dt + 'ij Qi{t)dX{t), Z(0) = 0, (2.19) 

where Qi{t) = w + ^4){t), (pit) = uj e°-'^du, and X{t) is a spectrally positive a-stahle Levy process 
with Laplace exponent — A°. 



Proof. By (|2.18|) . c^/n ~ °'L(cn) — > and c^/n ~ '^L{cn) — > 00, as n ^ 00. Then 
(C) holds. Without loss of generality, we consider G„(A) on A G [0,1] (see (I2.10p ). For A > 0, 
Gn{X) = -7A°[nL(c„)/<][L(c„/A)/L(cn)], and thus lim„_.oo Gn^A) = -7A°. For A = 0, the limit 
is trivial. Furthermore we have sup„ |G'„(A)| is bounded in A G (0, 1] as in the proof of Corollary 
12. li Note that GJ^(O) = and thus (D2) holds. It follows in a similar way that (E2) also holds 
and lim„_^oo -f^n(A) = — zzjA". By Theorem 12.21 the limit process Z{-) is described by ()2.15p with 
li{du) = (7a(a — l)/r(2 — a))u~^~"'du, u{du) = {'wa{a — l)/r(2 — a))u~^~°'du., g{-) = and 
Pi = [52 = 0. Note that Qi{t) > for t > 0. Define the process 

X{t)= / gi{sriuNoids,du)+ / / gi{s)-iul^o,<p{s)]{C)Ni{ds,du,dC). (2.20) 
Jo Jr+ Jo Jr+ Jo 

Then X{-) is a martingale. By Ito's formula, it is not hard to show that X{-) is a one-sided a-stable 
process with Laplace exponent —A". Thus we have (|2.19|) by (|2.2U|) and (|2.15|) . □ 



3 Asymptotic results for estimators 

In this section, we consider the statistical applications of our limit theorems as in [22^ lllj. For 
n £ N, suppose that a sequence of samples {(y„(A;), r?„(/c)), /c = 1,2, ••• ,n} is available. Then a 
natural estimator of the offspring mean m„ is given by 

n n 
"in = [^yn(fc - 1)] '^iVnik) -Vnik)). 
k=l k=l 

Using Theorem 12. H we can derive the following asymptotic result for m„. 
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Theorem 3.1 If the conditions of Theorem \2.1\ are fulEUed with F(X) being a unbounded function, 
and -R^(O) — > a as n ^ co for some a G M, then 

n{mn-mn) — > . , 3.1 

liY{t)dt 

where Y(t) and Y'{t) are dehned in Theorem \2.1[ 

Obviously, the above theorem apphes to the case of Coronary l2.1l Compared with the resuh of \22\ 
Corollary 3.1], this implies that the heavy-tailed stable distributions of offspring and immigration 
variables do not affect the rate of convergence of m„ in the critical GWI-process. 



We are also interested in the case when the conditions of Theorem 12.21 are satisfied. Then 
Yn{-)/n converges weakly to the deterministic function (p{t) = u f^e°'^ds, which implies only that 

n{rhn — TUn) — ^ by Theorem 13.11 Thus we further consider the applications of our fluctuation 
limit theorem, related to the CLSE of the offspring mean m„ based on only the information on 
{Unik)} as follows. For n,k > I let J^^ denote the cr-algebra generated by {yn{j)-,3 = 0, 1, • • • , k}. 
From (f2TT]) . 

nyn{k)\ri:_^] = mnVnik - I) + iOn. n>l. (3.2) 

If we assume that the immigration mean uj^ is known, then the CLSE m„ of m„, based on (|3.2p . 
is given by 

. _ Y2=iyn{k- l){yn{k) - UJn) , . 

If iOn is unknown, it is not hard to see that the joint CLSE {rhn-,0Jn) of {mn,ujn) is given by 



Ek=iiyn{k-i)-^] 



where 

1M = -y^yn{k), - Vy„(A:-l). (3.4) 



n ^ — ' n 

k=l k=l 



Using Theorem 12. 2^ we can derive the following asymptotic result for m^, and w^, which 
generalizes the result of |lll Theorem 3.1]. 

Theorem 3.2 If the conditions of Theorem \2.2\ are fulfilled with a; > 0, then 

where (/)(t) = lo £ e'^'ds, Z{t) is defined by ifZlg]) . M{t) = Z{t)-£ aZ{s)ds and it can be regarded 
as a deterministically-time-changed Levy process. Furthermore, 



— (m„ - m„ 




/ !^<i>(t)dM{t)-M{l) J^ct>{t)dt \ 

A/(l) ^Ht)dt-f^ 4>{t)dt 4,(t)dM(t) 



(3.6) 
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Corollary 3.1 Consider the CLSE rhn of rUn in the case of Corollary \2.3\ when condition (El) 
holds with u > 0. By the above theorem, we have 



n 



— (m„ - nin) — > U :-- 



where 1 < a < 2, c„ is given by li'J.lS]) . (j){t), Qi{t) and X{t) are given in Ii2.19]) . It is easy to see 
that U has a a-stable distribution and its Laplace transform equals 



E[e" 



exp 



f^cj)»{t)Ql{t)dt 



A > 0. 



Finally we turn to the case when the conditions of Corollary 12.21 are satisfied. In this case, it is 
possible to consider the CLSE estimates for the offspring and immigration variances -7r„ and 
Let Un{k) = Unik) — ninynik — 1) — ojn- Note that 



{k-l)+rn, n>l. 



(3.7) 



As in p6], if we suppose that nin and ujn are known, then the joint CLSE (vTn, '^n) of (7r„, r„), based 
on (j3.7p . is given by 



Efc=l(yn( 



^ul{k)/n-TTn 



(3.8) 



fc=l 



where y* is defined by (13. 3p . If ?ti„ and uon are unknown, we can use Unik) = ynik) — ThnVnik — 1) — 
Cjn instead of n„(A;) in (j3.8p and we get another joint CLSE denoted by {linifn)- Using Theorem 
2.21 again, we have the following asymptotic result for the above estimators, where the jumps of 
the fluctuation limit obviously play an important role. 



Theorem 3.3 If the conditions of Corollary \2.2\ i.e. (Dl), (El) and (a.1,2) are fulfilled with 
io > 0, then 



n(7r„, - TTn) 



' n ' n 



( mdj{t)-j(i) /; ^{t)dt \ 

j(i) <t>\t)dt-j^ <t>(t)dj(t) 
V /;02(iK-{/j<^(i)di)2 



where (j){t) = to e'^^ds, and J(t) is a martingale defined by 

ft f ft r rUs) 



J{t) =11 u''No(ds,du) + 

'o Jr+ jo 



u'^Ni{ds, du, d^). 



(3.9) 



(3.10) 



Jlere Nq and Ni are the compensated Poisson random measures given in i2.16\) . Moreover, i3.9\ 
still holds if TTn and fn are replaced by TTn and fn . 



Remark 3.1 A'^o and A^i are the Poisson random measures given in ()2.16p with intensities dsv{du) 
and dsfj,{du)dC, and v?v{du) + u'^fi{du) < oo. Let the limiting random vector in (j3.9p be 
denoted by (Ui, 1/2)'^ . It is not hard to see that if u^v(du) + u'^fi{du) < 00, then 



L- 



1 / /•! \ 2 ^00 

(^(/.(t)-y q){s)dsj dt J u^v{du) 

0(t) \y^{t) - J (j){s)ds] dt J uV(du) 
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where L = (fp'{t)dt — {f^ (j){t)dt)'^ . Otherwise we have E[C/^] = oo. Note that (f>{-) is not a const 
function. So if 7^ or /i 7^ (equivalently Nq or A^i is not degenerate), then Ui, and similarly 
U2, are not degenerate. 

We see that if the conditions of Corollary O are fulfilled with P(M+\{0}) > or /i(IR+\{0}) > 0, 
which means that the sequence of the offspring (or immigration) distributions fails to satisfy 
Lindeberg condition (b.2) (or (b.3)), then the resulting fluctuation limit Z{-) is a OU type process 
with positive jumps (see (j2.16p ). Thus, in this has the limit law Ui with normalizing factor 

n, and is not a consistent estimator. However, if we return to the case of [llj (see Example 
2.ip . which implies that Lindeberg conditions are satisfied and the resulting fluctuation limit Z{-) 



is a OU diffusion process without jumps (see ()2.17p ). then n(7r„ — 7r„) and Vn — rn ^ 0. In 
this case, to get the appropriate rates of convergence for 7r.„ and r„, we give the following theorem. 



Theorem 3.4 Consider the case of Example \2.1[ Let a^^n = E[(^„(l,l) — mn,)^] and b^^n 
E[(ryn(l) — i^n)^]- Suppose that (Dl), (El), (b.l) and the foUowing conditions hold with uj > 0: 

(c.l) na^^n — > 04 and b^^n — > 64 as n ^ 00 for some 04 > and 64 > 0, 



(c.2) nB (^„(1,1) -m„)^l|(g„(^i)_^^)2>^4 



as n — > 00 for all e > 0, 



(c.3) E {r]n{l) - UJn)'^'i-{(rjnil)-uJn)2>y^£} 



as n — > 00 for all e > 0. 



Let >P{t) = 00 £ e'^'ds, Q2{t) = 27r2(/.2(t) + (04 + 4^r)(/)(t) + (64 - r^) and V{t) = £ ^/g2{s) dW (s) , 
where W{t) is a one-dimensional Brownian motion. Let S — (^Jq (j^^it^dt — (Jjj (p(t^dt^'^^ (cij)2x2) 
where 

^11 = Jo ('^(*) ~ lo <P{s)dsY Q2{t)dt, (J22 = ( Jq (l)^{s)ds - (f>{t) Jq (f>{s)ds)'^Q2it)dt, 
(Ti2 = (J21 = Jo (</>(*) - Jo (A(s)ds) ( Jo 4>'^[s)ds - ^{t) J^ ^{s)ds)Q2{t)dt. 



Then we have 



ri^/^(vr„ -7r„,) 



( ^{t)dV(t)-V{l) <t>{t)dt \ 



\ 



V{1) cl>Ht)dt-f,] J^l Ht)dV(t) 



= AA(0,S). 

Furthermore, i3.11\) still holds if fcn and r„ are replaced by 7T"„ and f„, respectively. 



(3.11) 



Remark 3.2 It is easy to see that the above condition (c.l) implies that (b.2) and (b.3) in Example 
12.11 hold. So the conditions of our theorem is in the case of Example 12.11 If either of vr, 04 and 
64 — is not 0, then the limit normal law in ()3.1ip is not degenerate. 

Example 3.1 ( [111 Example 2.1]) The conditions of Theorem 13.41 are satisfied for the following 
examples with n(m„ — 1) — > a, n7r.„ a, na^^n cl and ?iE[(^„(l,l) — m„)^] — > a, as ?i ^ 00 
for some a > 0. (i) ^„(1, 1) has a Bernoulli distribution with mean 1 — an^^. (ii) the offspring 
distributions are geometric distributions with parameter pn = 1 — an~^, i.e. P(^n(l, 1) = i) = 



Pn(l - Pn 



A-l 



1,2,---. 
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4 Proof of main results 



Proof of Theorem \2.1\ For the proofs of Lemma 12.11 and 12.21 we can fohow the proof of Lemma 
2.31 or apply directly Corollary 1,2]. So we skip them. Now the limit functions R and F have 
representations (fTIl) and (123]) . Fix < A < M for any constant M > 0. Let \n = hn{l - e"^/''") 
and we have — > A. It follows from condition (A) that \Rn{\n) — Rn{^)\ ^ k{M)\Xn — A|, where 
k{M) > is a constant, and that lim„_»oo Rni^n) = -R(A). By condition (B) and the fact that Fn is 
a nondecr easing function on A G [0, M] for sufficiently large n, we have Fn ^ F locally uniformly. 
It implies that lim„^oo -^n(An) = -^(A). Let i?„(A) = i?n(A„) and Fn{X) = Fn{X). Note that the 
sequence {{Yn{ji),Yn{^)),l S N} is a Markov chain with state space En := {{i/bn,j/bn) : G 
N^} and the (discrete) generator An of t > 0} is given by 



= e-^'^'^^n exp { - 2;ia„(2;)e^i/^"-R„(zi)/n} exp { - /3„(2;)F„(zi + Z2)/n} 

= -e-<^'"> [Xianiz)e'''^'-Rn{zi) + Pn{z)Fnizi + Z2)] + o(l), 

where x e En, z = (21,2:2) > 0, q„(z) = (e'=i/*"5„(e-'^i/^") - l)"Mn (e^i/*"5„(e-^i/''")), and 
P^{z) = (/i„(e-(^i+^2)/bn)_i)-iln(/i„(e-(^i+^2)/bn))_ Ou thc othcr haud, Ict A bc thc iufinitesimal 
generator of {¥{■), ¥'{■)). For z » and x G^"^ 



^e-<^'^) = -e-<^'^) [xii2(zi) + F{zi + Z2)] • 

We need to prove that lim„_>oo sup^^^ \Ane~^^'^'^ — Ae~^^'^'^ \ = 0. The remaining proof is essentially 
the same as that of in |19l Theorem 2.1] or |2H Theorem 2.1] and so we omit it. □ 

Let us write / G C=k(M) if / is a bounded continuous function from M to M satisfying f{x) = o(a;^) 
when X ^ 0. Let F = [—1, cxd) and F„ = {(i — l)/cn : i G N}. Let fin be the distribution of ^• 
Then for sufficiently large n, ^„ is a probability measure on F supported by F„. 

Proof of LemmalM (sketch) Set 5n(A) = n2[e^^/^"(l - (m„ - l)A/c„) - 5n(e~^/'^")] and it 
follows from mean-value theorem that 

SniX) = Gn{X) + n^[mn - gMKe-^^"" - 1 + A/c„) 

W{1 - m„)(e-^/^" - 1 + X/cn) + n\mn - 1)(1 - e-^/^")A/c„, (4.1) 

where 1 - A/c„ <-!?„< e^'^/^". Under condition (D2), the sequence |G'„(A)| = ri^|s'^(l - A/c„) - 
iT^n\/cn is Uniformly bounded on each bounded interval [0, c] for c > and thus the sequence 
n'^\gn{'dn) — nT-nl/cn is also uniformly bounded. By (C), (Dl) and (D2), we have ^^(A) G{X) + 
^ajo, as n — > cx). To get ()2.12p . it is enough to consider the limit representation of Sn- Note that 



^y^-SniX) = -r? j (e-^^ - 1 + Xu)fin{du). 



We can use Venttsel's classical method (see [23] ) to prove it. More precisely, by modifying slightly 
the proofs of Proposition 2.1 and 3.1 in [20], we can show that there exist some constants (3i G K, 
(Ti > 0, and a c-finite measure ji defined as in (|2.12p such that 



(i) 



n 



(x('w) ~ u)fin{du) — > /3i as n — > 00; 
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(ii) 11^ / [u) ^,n{du) 2ai + / {u) ^{du) as n ^ oc; 

Jv Jo 

p poo 

(iii) lim / f{u)nn{du) = / f{u)fj.{du), for / G 
n-^oo _/p 7o 

Note that e"^^ - 1 + Ax(x) - ^X'^x^ix) G C^(r) as a function of x G T for fixed A > 0. The 
above results imply that the limit function of Sn has a Levy-Khintchine type representation. Let 
(Ti = ai + ^ajo and let /?i = /3i + J^{u — x{u))f^{du). Then we have (j2.12p . But we still need to 
verify ui > 0. It follows from (a.l), (C) and (Dl) that 

— / x{u)l^n{du) = — / {x{u) - u)f,nidu) + (4.2) 

which tends to 070 as n — > 00. Let E be the set of e > for which fJ,{\u\ = e) = 0. By ()4.2p . (ii) 
and (iii), we obtain 



lim lim [ ( x\u) + ^]iin{du) = 2ai. 



The support of iin is r„ and for large enough n, x^{u) + (x(^) /cn) > if n G r„. Thus ai > O.D 

Let f ,„ = {ijcn : « G N} and let Vn be the distribution of Then Vn is a probability measure 

on [0,00) supported by r„. 

Lemma 4.1 Under conditions (C), (El) and (E2), <l2.i3|) holds. As n ^ 00, we also have 

(i) n / {x{u) - u)vn{du) ^ (32 - / (u - x(^^))j^(rf^^); 

(ii) n / {u)^n{du) ^ 2a2 + u:-^Q + / x^{u)i^{du); 

Jo Jo 

roo roo 

(iii) limn/ f{u)vn{du) = / f{u)v{du), for / G C^(]R+). 

Jo 7o 

Proof. This lemma is proved with the same method as Lemma 12.31 But we need to prove that 
2(72 + uJ^o ^ "^^70 • Let dn = x{u)vn{du) and let E be the set of e > for which v{u = e) = 0. 
By (C), (El), (i) and (ii), it is not hard to show that 



lim lim n [ (x^(^) ~ ^ )'^n{du) = 2a2, 
lim lim n / ix{u) - a„)^f„((iti) = 2(72 + ^7o - t^^To- 



For large enough n, x^('w) ~ (x(^)/cn) > if n G Tn- Then we are finished. 
Lemma 4.2 Under the conditions of Theorem \2.2\ we have for t > 0, 

lim sup — E 

n—*oo ri 



sup y„(s) < \a\<P{t) +ujt, (4.3) 

0<s<t -I 
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limsupE sup |Z„(s)| < M(t) [l + (|a| + l)t exp{(|a| + l)t}] , (4.4) 

where = /g e^\''\+'^> duds , M{t) = 2K{^{t) + t) + A^yK{'^{t)+t), and K is a positive 
constant defined as in H.9]) . 



Proof. Note that (12. ip can be rewritten into the following form: 

I yn{k-l) I 

yn{l)=Yl E iUKj)-l)+Y.'^n{k). 
k=l j=l k=l 



(4.5) 



Let Uk,j) = iUkJ) - l)/cn, wn{k) = ZT=i ixiUkJ)) - B[x{Uk,j))]), and Wn{l) = 
Ylk=i'^rL{k). let denote the u-algebra generated by {{wn{j) , yn{j)) , j = 0, 1,-- - ,k}. Since 
E[wn{k)\T]^_^] = 0, VFn(["'i]) is a square integrable martingale, and the quadratic variation is 
Yl^k=i'^nik)- On the other hand, it follows from conditions (Dl) and (El) that 



1 f[ 



nt]/n /•[ns]/n 



m. 



l^^^'^duds, 



(4.6) 



which tends to lo Jq e°''^duds, as n ^ oo. Then applying Doob's inequality to martingale terms 
in ()4.5p . we have for sufficiently large n. 



E 



sup Yn{s) 
0<s<t 



+ 2c„E2 



[nt] y„{k-l) 

<CnB[Y^ {Uk,j)-x{Uk,j))) 
k=i j=i 

[nt] 
k=l 



+ n^c„$(t) / x{u)fj.nidu) 



+ nujnt- 



By (i), (ii), (C), (Dl) and (El), we obtain ncn x{u) Hn{du) — > a and then (j4.3p holds. By (|4.5p . 
the sequence Zn{-) are given by 

[nt] [nt] y„{k-l) [nt] 

Zn{t) = }_^{mn-l)Zn{—-)+}_^ + • (4.7) 

k=l k=l j=l k=l 

Let f]n{k) = rjn{k)/cn- By Doob's inequality, it is not hard to see that for sufficiently large n, 



E 



sup |Z„(s)| <n|m„-l|/ B[\Zn{s)]ds + 2n'^t) f {u - x{u))^Jin{du) 

0<s<t -I 



r 

oo 



2 n^$(t)varx(en(l,l)) + 2nt {u - xiu))iyn{du 



+2(ntvarx(??n(l)) I • 



(4.8) 

By Gronwall's inequality and standard stopping argument, (i), (ii), (C), (Dl) and Lemma 14.11 
implies 



E[\Zn{t)\] < 2{Kmt) +t) + 2{K{^{t) + t)) 2 } exp{(|a| + l)t}. 



(4.9) 



14 



where K = sup„ (n^ /r(''^ ~ x('w) + {u)) fin{du) + n f^iu — x{u) + x^{'u))i^n{du)) . By the above 
inequahty and (|4.8p . we obtain (|4.4p . □ 

Lemma 4.3 Let 0n(i) = E[y„(t)]/n for t > 0. Under the conditions of Theoreni \2.2\ the sequence 
{Zn{-),M-)) is tight in D([0,oo),M x M+). 

Proof. By Lemma|121 C{t) := 1 + hmsup„^oo (iE[supo<s<t y„(s)] + E[supo<^<( |^n(s)|]) is a 
locally bounded function of t > 0. Then Zn{t) is a tight sequence of random variables for every 
t > 0. Now let {r„} be a sequence of stopping times bounded by T and let {Sn} be a sequence of 
positive constants such that (5„ ^ as n — > 0. By Doob's Optional Sampling Theorem, we obtain 
as in the calculations in (j4.8|) that for sufficiently large n, 

B[\Zn{Tn + 5n)- ZniTn)\] 

< 2K / " -E[yn{[nTn] + [ns])] ds + (|a| + 1) / " E[Z„( ^^"""^ ^ ^^^^ )] 
Jo n Jo n 

lnSn] + l 1 - 

+ (k " -B[yn{[nTn] + [ns])] ds) ' + 2K(5n + -) + + - 

\ jQ Tl / \ 71/ \ Tl 



<(2K + |a| + l)y C{T + s)ds+{K j C{T + s)ds^ 

+2K(5n + -) + (K{5n + - 
\ nJ \ n 



Then Zn{-) is tight in L'([0, oo),M) by the criterion of Aldous [Tj. It is easy to see that (l)n{t) 
converges to (t){t) := uj j^e'^'^ds in distribution on Z)([0, oo), M+). By Jocod and Schiryaev [TH 
Corollary 3.33, P.317], (Z„(-), </>«(•)) is tight in Z)([0,cx)),M x IR+). □ 

Let Z{-) be any limit point of Zn{-). Without loss of generality, by Skorokhod's theorem, we 
can assume that on some Skorokhod's space {Q,J^,J^t,P), {Zn{-), 4>n{-)) {^{')i4>{')) the 
topology of L>([0, oo), M x M+). 

Lemma 4.4 For any fixed A E M, 

L{t) = e^^^W _ eiAZ(o) _ /"* e^^z{s) ^^^^^^^ ^^^^^ ^^^^ (4_-^q) 

Jo 

is a complex-valued local J^t-^8,rtingale. Here = — 1 and 

A{xi,X2, A) = iaXxi + (a7oA^/2 - G{-iX))x2 + 7o(u;^ - a;)A^/2 - H{-iX), 

where G and H are defined hy Ii2.12\) and \2.13\) . respectively. 

Proof. Define the stopping times 

/ = inf{t > : \Z{t)\ > 6 or \Z{t-)\ > b}, 
tI = inf{t > : \Zn{t)\ > b or |Z„(t-)| > b}. 

Let Z\t) = ZitAr''), Z^{t) = A r^), and analogously /(t), (/>^(t). It follows from [121 

Proposition 2.11, P.305] that for all but countably many b, 



h a.s. h 



r^nR and (Z^(.), <^:;(.)) ^ (Z''(.), <^^(-)) 
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in the topology of D{[0, co),M. x M+). Define T^{t) = M and r*(i) = M. We claim that 

r^(-)^r''(-) in C([0,oo),M+), as n ^ oo. (4.11) 

In fact, since < T^(i + e) — T^(t) < e for any t > 0, the criterion of Aldous yields tightness for 
{<(•), n > 1}. On the other hand, {^n(|) : A; > 1} is a time-inhomogeneous Markov chain. For 
fixed A G M, 

fc=0 

is a complex-valued martingale. ()1.2p implies that 

L„([nt]) = e*^^"W - e*^^"W - / " e'^^-^'^n[An[Zn{s),Ms)A) - l\ds, (4.12) 

Jo 

where An(xi,X2,A) = e-^V^n (n(m„-i)x2+<^„)(g-a/c„^^^giA/c„)^c„xi+nx2^^^giA/c„)^ Pqj. simplicity, 
we denote L„([nt]) by Ln{t). Then L^(t) := Ln{t A r^) is also a complex-valued martingale. It 
follows from the proof of Lemma 12.31 and Lemma 14.11 that 

n(e-^^/^"g„(e*^/^") - 1) ^ and n^hn{e'^/^") - 1) ^ ioj-fl X, (4.13) 
as re ^ oo. Then we have for sufficiently large re, 

/>oo 

+(c„xi + rex2)/i,n(A) + / (e*^"* - 1 - iXu)un{du) 

Jo 

-\{hn{e'^/'")-l)'+l2,nW, (4.14) 



where 

pjA/c„^_('giA/c 



/i.(A) = g(-i)^-^[--^^"^"^7"^-^^ 

/..(A) ^ f;(-i)-^M^:^^!:l^. 



Note that n2|/i^„(A)| < |re(e*^/^"c/„(e*^/'=") - 1)^ ^ and re|/2,„(A)| ^ 0. By (i)-(iii), LemmaSH 
(j4.13p and (j4.14p . it is not hard to show that n{An{xi, 2:2, A) — 1) — > A{xi, X2, A) locally uniformly 
on {xi,X2) G M X for fixed A. As in Ethier and Kurtz O Problem 26, P153], we obtain that 



* e*^^"(^)re[A (Z^(5), 0^(s), A) -l]ds^ f e'^^' '^'^ A{Z\s) , <p\s), \)ds 



in the topology of C([0,oo),C). Let L^{t) = L{tAT^). Note that [ret]/re t in C([0, 00), 1R+). By 
(liTTT) . [3 Problem 13, P.151] and [121 Proposition 1.23, p.293], we have 



L^(t) ^ L''(i) in D{[0, 00), C), as re ^ cx). (4.15) 

Then for almost all t > 0, L^^i^) -^*(^) C!. Fix arbitrary T > 0. For any t < T, 

I /o""^*^ e*'*'^"'-*-* Z^(s)(is| < 6r, where the bound holds uniformly in n. Then for almost t < T, 
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-^n(^) as n — > oo. Since Lf'{t) is right continuous and bounded for t <T, we have 

is a martingale. Note that r'' — > oo as 6 ^ co, L(t) is a local martingale. □ 

It follows from ()4.10p and \T2[ Theorem 2.42] that Z{-) is a semimartingale and it admits the 
canonical representation 

Z{t) = Z(0) + + / [132 + (3i(j){s) + aZ{s)) ds + j j uJ{ds, du), (4.16) 

Jo Jo Jo 

where Z{0) = 0, Z'^{t) is a continuous local martingales with quadratic covariation process Jq g{s)ds 
with q{s) = (2(Ji — a^Q)(p{s) + 2a2 + uj{l — 0^)70, and J{dt, dz) is an integer- valued random measure 
on (0, 00) X R_|_ with compensator J{dt, du) = (j){t)dt^{du) +dtv{du), where J{dt, dz) = J{dt, dz) — 
J{dt, du). 

Lemma 4.5 Suppose that the conditions of Theorem \2.2\ are satisfied. Then the cadlag process 
Z{-) is a weak sohition of i2.15\) . 

Proof. Define the measure p{du,d(^) = ^{du)i{dC,) + u{du)6Q{dC,), where i{dC,) is the Lebesgue 
measure on (0, 00) and 5o{dC,) is the Dirac measure at = 0. By Ikeda and Watanabe [ini P-84 
and P.93], there exists a standard extension of ($7, Tt, P) supporting a one-dimensional Brownian 
motion and a Poisson random measure N{dt, du, dC) on (0, 00) x with intensity dsp{du, d() such 
that dZ^it) = y^~Q(t)dB{t), and 

J{{0,t]xE) = f [ lE{e{s,u,C))N{ds,du,dC), (4.17) 





for any E G 53(R+), where d{s,u,C) = ul[o,<^(s)] (C)- Set No{ds,du) = N{ds,du,{0}) and set 
Ni{ds,du,dC) = A^(c?s, dw, (iC)l(o,oo)xiR4.x(o,oo)- Then we see that Z{-) is a solution of ()2.15p . □ 

Proof of Theorem \2.2\ By [IQi P.231], the Lipschitz conditions of the equation ()2.15p imply its 
pathwise uniqueness of solutions. Thus Theorem 12.21 follows from Lemma 14.31 and 14.51 □ 

Proof of Theorem\3j\ By Theorem O {Yn{-) /bn,Y^{-)/bn) converges weakly to {¥{■),¥'{■)) 
on L'([0, 00), M^), and (¥{■), Y'{-)^ is stochastically continuous. Note that 



mm — mr, 



Yn{l)/hn - n{mn - 1) Yn{t)/bn dt - Y^il)/bn 
fiYn{t)/bndt 

Then we have (13. ip by the continuous mapping theorem. Since F is not bounded, the immigration 
process Y'(-) is neither a compound Poisson process or a zero process. This implies that P(Y{t) = 
for alH G [0, 1]) =0. □ 

Lemma 4.6 Define Un{k) = yn{k) — mnUnik — 1) — w„. Tiien we have 

1 " 

y"ul(k)^0, asn^oo. (4.18) 

ncn 

" k=\ 

Proof. It follows from and ([T^ that 

'"n(^) = iCn{k,j) - m„) + {r]n{k) - Un)- 

i=i 
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Recall that Cn{k,j) and f}„(/c) defined in the proof of Lemma[12J Note that Cn{k, j)—x{S,n{k, j)) > 
and fin{k) — xiVnik)) > 0. Then we have 

ncn n L a Jq \ n / 

+ 6c„f / {x{u) - u)Vn{du)\ \ —^{fln{k) - Xifjnik))) 



k=l 



6c 



(4.19) 



k=l 



where 



^ k=l j=l 

n y„{k-l) 

E U(ln(^,j))-E[x(|„(A:,j))]) 



'l,n 



'2,n 



fc=l i=l 



We obtain that E[/i,„] < y^<^{l)K and E[/2,n] < ^$(1) [n^ var x(|n(l, 1))] as in the calcula- 



tions in ()4.8p . Condition (C) implies that as n ^ oo for i = 1,2. From (C), Remark 



and (a.l) in the proof of Lemma |2.3[ the third term in (|4.19p converges in probability to as 
n — > oo. As in the above proof, we also have that the last three terms converge in probability to 
0. Thus implies (ITO]) . □ 

Proof of Theorem \3.S\ First consider the equation ()2.15p . We obtain as in the calculation in 
(giH]) and (filH that for < s < t, 

E[\Z{t) - Z{s)\] < |/32|(t-s) + |/9i| / 4>{u)du + \a\ j $(u)el"l"dn 

J S J S 



+ 2/ 4){u)du j {uAu^)fi{du) + 2{t-s) {uAu^)u{du) 



+ 



u)du I {uAu )i^i{du)] + I / Q{u)du 



(4.20) 



where ^>(-) is some non-decreasing continuous function. Then Z{-) is stochastically continuous. 
Let D{Z) := {t > : Y>{Z{t) = Z{t-)] = 1} and thus D{Z) = (0,00). From we obtain 



n 



{rhn - nin) 



_J_ v^n 

nc 



Y2=iyn{k - l)Un{k) D{ 



n 



^Trk=iyl{k-i) Q{ny 

Rewrite D{n) as D{n) = Di{n) + ^ Y^^j=2 Dj{n) - Di{n), where 

Di{n) = —j2nyn{k-l)]un{k), I)2(n)-''^^""'"^ 



(4.21) 



fc=i 



2m- 



n JO 



Zn{s)ds, 



Let Mn{t) = Y}k=i'^n{k) / Cn for t > 0. The functional ■ ^([0,oo) 



2m„ 



A(n) 



^E-^(^)- 



2ncnmn 



(4.22) 



k=l 



is defined by 



^n{x) = UJn [\x{l) - X{t))m^;:'^-'dt - 



(4.23) 
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Then Di{n) can be rewritten as 



^ n—l n 

Di{n) = V Un{k) = ^'„(M„ 



j=i k=j+i 

If Xn ^ X in the topology of D{[0, C!o),M), it is easy to see that |^'n(x„) — ^(x)| — > 0, where 

-1 

at , 



^{x) =cj / (x(l) - x(t))e''*dt. 

Note that M„(t) = - /(|"*^^"n(m„ - l)Z„(s)ds by (gZ]). It follows from Theorem 

that Mn{t) converges weakly to M{t) := Z{t) - jlaZ{s)ds on i:>([0, oo), M). By Remark 
Yn{-)/n converges weakly to on D([0, oo), M4.), and cj){-) is a deterministic continuous func- 
tion. Thus (y„(-)/n,Z„(-),M„(-)) converges weakly to , Z {■) , M {■)) on L>([0,oo),M+ x M^). 
By Theorem 7.8, P. 131 and Problem 26, P. 153] and the continuous mapping theorem, we have 

Di{n) 4 ^'(M) = /oV(s)dM(s), 4 - /^^ aZ(s)ds, Z^sW ^ \Z'^{1) and Q„ ^ /oV(s)ds, 

as n — > 00. Then it follows from (C) and Lemma 14.61 that ^ Y^j=2 ^ji^) ~ -^4("') ~^ • Hence we 
obtain (13. 5p . In a similar way, we also have ()3.6p . □ 

Recall that tin(A;) = yn{k) — mnUnik — 1) — uJn- Let Vn{k) = u^i{k) — 'KnUnik — 1) — r„ and let 
= Y}k=iM^)- By (EZI), Vn{-) is a martingale. 

Proof of Theorem VJ.'J\ Under the conditions of Corollarv 12.21 Zn{-) is defined by (|2.14|) with 
Cn = y/n, and then Mn{t) = Yl'h^i Un{k) / y/n. By Corollary 12.21 and the proof of Theorem 13.21 we 
have that (Z„(-), M„(-)) converges weakly to (Z(-),M(-)) on D([0, cx)), M^)^ where Z(-) IS given 
by ()2.16p and M{t) = Z{t) — aZ{s)ds. Note that M„(-) is a square integrable martingale 
and B[M^{t)] = (vr„, E[y„(A; - 1)] + [nt]r„)/n. Then for i > and sufficiently large n, 

B[M^{t)] < /g (j){s)ds + P(M+)t + 1. Thus by Kurtz and Protter [El Theorem 2.7], 

M„(t), Zn{s-)dMn{s)j ^ (^Z(t), M(t), ^* Z(s-)a!M(s)) (4.24) 

in distribution on D{[0, cx)),M^). On the other hand, let Vn{t) := Vn{t)/n. 

Vn{t) = Zl{t) + n(l - ml) / Zl{s)ds - 2m„ / Zn{s-)dMn{s) 

Jo Jo 

r-[nt]/n 

-niTn Yn{s)/nds - [nt]rn/n. 



Still note that y„(-)/n converges weakly to (/>(•) on D([0, cxd), M+), and is a deterministic 
continuous function. By (14.240 and the continuous mapping theorem, (Z„(-), V^(-)) converges 
weakly to (Z(-), J(-)) on L>([0, 00), M^), where J{t) = Z^{t) - 2a Z^{s)ds - 2 Z(s)dM(s) - 
/i(M+) Jq (j){s)ds — D(R^)t. By Ito's formula, J{t) has also the form (|3.10p . Vn{t) is also a finite 
variation process. Denote its finite variation by /q |(il^(s)|. Then for t > and sufficiently large 
n, B[ \dVn{s)\] < 2/2(M+) £ (/>(s)ds + 2P(M+)t + 1. By [l5] again, £ y„(s-)/ndK(s) converges 
weakly to £ (l){s)dJ{s) on D([0, 00), M). We see that 

, fo Ynis-)/ndVnis) - Vnjl) li Yr,{s)/nds 

nlTTn — 7^n) = — T 1 , (4.25) 

Jo\Y4s)/n-J^'Y^{s)/ndsyds 
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and f„, — r„ = t4(l) — (^Tn — vr„) J"J y„(s)ds. Note that J{t) and (j){s)dJ{s) are stochastically 
continuous. By the continuous mapping theorem, we have (j3.9p . We write 

'U„(fc) - = -[{rhn - mn)yn{k - l)f - (a)„ - lOnf + 2(m„ - mn)yn{k - l)un{k) 

+2{Qjn - iOn)Un{k) - 2(m„ - ?n„)(a)„ - uJn)yn{k - 1). 

As in the proof of (j4.25p . also by Theorem 13.21 we have (j3.9p holds for 7f„ and f„. □ 

By the proof of Theorem 13.31 we see that Vn{t) := Vn{t)/n converges weakly to J(-) on 
L'([0, oo), M), where J{t) is defined by (|3.1Up . However when we turn to the case of Example 
12.11 J{t) is degenerate to 0. Then in this case, we need the following lemma. 



Lemma 4.7 Let Vn{t) = Vn{t)/\/n. Under the conditions of Theorem \3.4l Vn{-) converges in 
distribution on Z)([0, oo), M) to the process V{-), which is defined by Vit) = ^/g^dWis), 
where Q2{t) = (f)^ {t) + (04 + 47rr)0(t) + (64 - r'^) and W{-) is a one-dimensional Brownian 
motion. 



Proof. Under the above conditions, we see that y„(-)/n converges weakly to </>(•) on -D([0, 00), M+) 
by Remark 12.31 (ii). Then for any t > 0, 

^ [ni] ^ [nt] 

- J2 nvlik)\^Li] = - E ((«4,n + 47r„r„ - 3nl)yn{k - 1) + 27r2y„(fe - 1)2 + 54,^ - rl) , 
n ^ — ' n ^ — ' V / 

fc=i fc=i 

which converges in probability to Jq g2{s)ds. Now by the martingale central limit theorem, it 
suffices to prove that, for any e > and t > 0, 

^ [nt] 

-E^[^"(^)l{K(fe)|>v^e}l-^-i] ^0, asn^oo. (4.26) 

k=l 

We have that Vnik) = An,k + -Bn.fc + Cn,k + Dn,k, where 

yn{k~l) y„{k-l) 
An,k = ^ [{Cn{k,i) - m,nf - VrJ, BnM = 2 ^ {in{k,i) - mn){r]n{k) - UJn), 

1=1 1=1 

yn{k-i) 

Cn,k = (rinik) - UJnf - rn, Dn,k = 2 ^ {Cn{k , i) - m.n) {Cn{k , j) - ■ 

i<j 



Note that for any pair of random variables X and Y, 'E[{X +Y)'^l^^x+Y\>e}] — '^y^[-^'^^{\x\>£/2}] + 
E[y2i^l^l^^^2|]^. Thus it suffices to show that (|4.26p with Vn{k) replaced by An^k^ Bn.k^ Cn,k, 
and Dn,k- Let ^'^{k,i) = {£,n{k,i) — mn)^ — vr„. As in the proof of [HI Theorem 2.2], we obtain 

^ [nt] 

~'^^\^n(^)'^{]Ar,(k)]>V^e}\H-l\ 
k=l 



[nt] [nt] 

< nE[(C;(l, l))'l{|5;(i,i)|>VH./2}] Yl y-^^ - 1)/^' + ^"(«4,n - vr^) V2 ^ yl{k - 1)/ 

k=l k=l 

[nt] 3 

+V2n{a4,n - 7r^)^e"^ ^y^ik- l)/nt. 

k=l 
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For large enough n, 7r„ < and 

E[(Cn(l, l))^l{|5;(l,l)|>^e}] < E[(^„(l,l) - m„)^l|(g^(^i)_^^)2>^e/2}] 

+2(7r„a4,„ + 2'K^^|^fne. 

Then, by conditions (c.1,2), (j4.26|) holds with f„,(A:) replaced by A^^k- Also by condition (c.3), 
()4.26p holds for Cn^k- Let (,nik,i) = S,n{k,i) — rUn- For Dn,k, we note that 

j/„{fc-l) y„{k-l)-l 

= vr„ Yl (j-l)C(^,j)+vrn Yl {yn{k-l)-iKnik,i) 

j=2 i=l 
yn(fc-l) y„{k-l) 

+ E C(A:,i)C;(A:,j)+2 J] {Uk,l)fUk,i)Uk,j) 

i<j l<i<j 
y„{k-l) y„(k-l) 

+2 Uk,l){Uk,i)?Uk,j)+2 Y Uk,l)Uk,iMn{k,j)? 

l<i<j l<i<j 

yn(k-i) 

+6 Y UkJ)UkJ)Uk,j)Uk,p) + yn{k-l){ynik-l)- 1)7^^/2. 

l<i<j<p 

Then it follows from the above equality that 

E[<,|^_i] < Waiy^ik - 1) + 416a4,nvr2y3(fc _ i) + 772<y4(/, _ i). (4.27) 

Thus, for any t > 0, ^ ES E[D2(A:)l^p^(,)|>^,}|^^ J < ^ EHi E[i^^(A;)|^,"_i], which 
converges in probability to by (14.270 . In a similar way, we can also prove that ()4.26p holds with 
Vn{k) replaced by i?n,fc- D 

Proof of Theorem \3.4\ It is not hard to see that for any t >0, 

[nt] [nt] [ns] lnu]_ 

-sY^iyUk-l)] = ^Ii±^ / " / " ^-M / " rr^inCl,^ duds 

" ^ "^n JO Jo Jo 

[nt] [ns] 



Jo 



which converges to fQ(jP'{s)ds as n ^ 00. Also by (j4.6p and the proof of Lemma 14.31 we see 
that Vn{-) is a square integrable martingale and sup„E[y^(t)] < 00 for t > 0. By [151 The- 
orem 2.7], [Yn{t)/n,Vn{t), jlYn{s-)/ndVn{s)) [(t){t),V{t)Jl(t){s)dV{s)) in distribution on 
Z)([0,oo),M+ X m?). As in the proof of Theorem [331 we have (|rm . □ 
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